Crossover from negative to positive shear rate dependence in granular friction 
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We conduct an experiment on the frictional properties of granular matter over a wide range of 
shear rate that covers both the quasistatic and the inertial regimes. We show that the friction 
coefficient exhibits negative shear-rate dependence in the quasistatic regime, whereas the shear-rate 
dependence is positive in the inertial regime. This crossover from negative to positive shear-rate 
dependence occurs at a critical inertial number. This is explained in terms of the competition 
between two physical processes, namely frictional healing and anelasticity. We also find that the 
result does not depend on the shape of the grains and that the behavior in the inertial regime is 
quantitatively the same as that in numerical simulations. 
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The rheological properties of granular matter are im- 
portant for analyzing various phenomena in geosciences, 
such as landslides, debris flows, and crater formation. 
These properties may also be relevant to friction in faults 
because a fault contains fine rock powders that have been 
ground up by past fault motions. Among the rheological 
properties of granular matter, the shear-rate dependence 
of the friction coefficient determines the stability of de- 
formation and therefore has been paid major attention. 
During the last decade, considerable progresses have been 
made on this subject; it is recognized that a nondimen- 
sional number plays a central role in describingthe shear- 
rate dependence of the friction coefficient . This is 
referred to as the inertial number ^ = 7 a/ m/ Pd, 
where 7 is the shear rate, P is the pressure, m is the 
mass of the grains, and d is the diameter of the grains. 
As m/ Pd is the velocity relaxation time in granular 
matter 4], the inertial number may be regarded as the 
index for the inertial effects inside the flow and there- 
fore one can classify the flow regimes according to the 
value of the inertial number. The flow is fast and col- 
lisional at larger inertial numbers (/ > 10^^), which is 
referred to as the inertial regime. At lower inertial num- 
bers (/ < 10""'^), which is conventionally referred to as 
the quasistatic regime, the grains are closely packed and 
slowly sheared. In the inertial regime, a constitutive law 
obtained in numerical simulation quantitatively explains 
the experimental observation on an inclined plane flow 
0. However, we have not reached the consensus about 
a constitutive law for the quasistatic regime [^-Q. Al- 
though some numerical studies 0, 0-11 1 are conducted 
down to / ^ 10~^ to address a constitutive law for the 
quasistatic regime, they are yet to be verified in experi- 
ments. In addition, we are not aware of any constitutive 
laws that can lead to the exponential veloci ty p rofile that 
is ubiquitously observed in heap flows [13, [l3| . 

To clarify the nature of granular friction in the qua- 
sistatic regime and its connection to the constitutive law 
in the inertial regime, one must conduct the rheologi- 



cal measurements over a wide range of shear rates that 
cover the both regimes. To this end, we used a com- 
mercial rheometer (AR-2000ex, TA Instruments), with 
which we can control the sliding velocity over a wide 
range (1 x 10"^ — 3 m/s). The normal force is also control- 
lable. We chose a relatively low pressure (20 — 50 kPa) 
to exclude the effects of frictional heat In addition, 
the temperature beneath the lower plate was set to 25° C 
and was kept constant with Peltier Plate. To study the 
effects of the grain shape on the rheological properties, 
we prepared two kinds of grains: soda-lime glass beads 
(spherical with a mean diameter oi d = 270 fim) and 
chromite sand (angular shaped with a mean diameter of 
d = 286 fim). The grains were packed into the annular 
channel of the sample holder [15| as shown in FIG. [1] 
a. The cylindrical side wall, which was fixed to the lower 
plate, is made of transparent fused silica so that we could 
optically observe the internal structure of sheared granu- 
lar matter |16| . The images of each grain were captured 
by a high-speed video camera (Photron Fastcam APX 
RS). 

To observe the internal velocity profile, we used parti- 
cle tracking velocimetry to obtain the instantaneous ve- 
locity of each grain, which was then averaged out with 
respect to the flow direction. The velocity proflles are 
shown in FIG.[T]b. Note that the velocity proflles that are 
normalized by the sliding velocity V collapsed to a single 
curve. As the velocity profiles were somewhat nonlinear 
as observed in we defined the effective shear rate 
7 = V/Ws using the effective flow width Wg, which was 
deflned as the depth at which the normalized flow velocity 
decreased to 1/10. We estimated Wg — 5d. Throughout 
this study, the inertial number is deflned in terms of the 
effective shear rate. 

The shear stress a and the layer thickness H were mea- 
sured at each steady state of the sliding velocity V under 
a constant normal stress P. Fig. [l]c shows an example 
of the normal-stress dependence of the steady-state shear 
stress. The intercept of the best-flt line is quite close to 
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FIG. 1. (a) Schematic diagram of the apparatus. A gran- 
ular sample fills the annular channel, which is attached to a 
rheometer. Shear is applied to the sample by rotating the up- 
per plate, whereas the lower plate is fixed. To avoid slipping 
on the surfaces, a monolayer of the glass beads was glued to 
both the upper and the lower plates, (b) The velocity profiles 
normalized by the applied sliding velocities, V, ranging from 
10"'' to f m/s. (c) The normal-stress dependence of the shear 
stress. 
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FIG. 2. The friction coefficient plotted against the inertial 
number. Symbols show the different normal stress (light blue 
circles: glass beads 20 kPa, blue circles/triangles/squares: 
glass beads 30 kPa, purple triangles: glass beads 50 kPa, 
green triangles: chromite sand 30 kPa). 



the origin and therefore the resistance from the side walls 
may be neghgible. We then define the friction coefficient 
/J, as a /P. In Fig. [51 the friction coefficient is shown as 
a function of the inertial number. 

At lower inertial numbers (/ < 10^^), the friction co- 
efficient is weakly dependent of the sliding velocity. Note 
that the dependence is negative and logarithmic. The ex- 



tent of weakening is characterized by the slope d^/d\nl. 
For example, the slope is approximately —0.003 for spher- 
ical glass beads at 30 kPa. The negative dependence on 
the sliding velocity is indeed well known in earthquake 
physics and in tribology [l^ and the extent of weak- 
ening in the present study is comparable to those typi- 
cally observed for glass or rocks, although the slope may 
be sensitive to the experimental details, such as humidity. 

Although the negative slope is the common feature of 
the present data, the absolute value of the friction co- 
efficient differs significantly from sample to sample. An 
apparent ingredient that affects the absolute value is the 
grain shape [2^ [2l| . The friction coefficient of the angu- 
lar grains (/Lt ~ 0.55) is significantly higher than that of 
the spherical grains {fj, ~ 0.4). It is also noteworthy that 
the spherical glass beads may exhibit anomalously low 
friction (/i ~ 0.3) despite the same experimental condi- 
tions (e.g., the blue squares in Fig. 2a). This effect may 
be due to the structural ordering induced by the shear 
[2^ . although that hypothesis has not been verified. 

At higher inertial numbers (/ > 10^^), the friction co- 
efficient substantially increases and the data indicate the 
characteristic inertial number above which the friction 
coefficient increases. We refer to this as the crossover 
inertial number, and it is denoted by Ic- Although Ic 
for each sample varies slightly, its value is apparently 
on the order of 10^^. Thus, for simplicity, we chose 
Ic — 0.032 assuming that Ic is common to all of the 
data. We then defined the amount of strengthening as 
Afi{I) = /i(/) — /i(/c) and found that A/i(/) collapses in 
the high / regime (/ > Ic) (as shown in Fig. [3^). This 
collapse obeys 



A^(/) cx ci/. 



(1) 



with ci ~ 0.6. This strengthening behavior for I > Ic 
is accompanied by dilation. To compare the data with 
different layer thicknesses, it is convenient to define the 
nondimensional dilation as AH* = {H{I) — H{Ic))/Ws- 
As shown in Fig. [3]b, AH* also collapses in the high / 
regime. 

AH* cx C2I, (2) 

where C2 — 0.2. From Eqs. H]) and it follows that 

Afi = C3AH*, (3) 

where C3 = ci/c2 — 3. Figure [5] shows that Eq. ([3]) 
explains the experimental data well. In addition, we no- 
ticed that Eqs. (P), an d (Bt quantitatively explain 
the inclined-plane flow data [23|. Note that the propor- 
tional coefficients (ci, C2, and C3) are common to these 
different experimental data. We thus conclude that these 
relations are independent of the details of granular mat- 
ter such as the grain shape and the coefficient of restitu- 
tion. 
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FIG. 3. (a) The increase of the friction coefficient, A/i(J)(= 
fj,{I) — fi{Ic)), where Ic{= 0.032) is the crossover inertial num- 
ber. The inset is the linear plot, (b) The normalized thick- 
ness change AH* (I). Symbols are the same as those in the 
previous figure. The inset is the linear plot. Results of the 
DEM simulation with an inter-granular friction coefficient of 
0.6 and a restitution coefficient of 0.8 are also shown (black 
diamonds) . 
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FIG. 4. Friction coefficient difference Ajj, versus normalized 
thickness change AH*. The dashed straight line indicates a 
slope of 3. Symbols are the same as those in the previous 
figures. 



Relations (P), 0, and Q, which hold for the high / 
regime, are quantitatively the same as those seen in nu- 
merical simulations P, 0, [1, 0-[li|L The data obtained 
from the numerical simulations [9| are also shown in 
FIGS. [3] and m which show that the simulations quanti- 
tatively agree with the experiments in the high / regime. 
It is believed that the velocity strengthening seen in nu- 
merical simulations is caused by the anelasticity of the 
grains, which is usually modeled by a "dashpot" 
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The quantitative agreement between the simulation and 
the experiment in the high / regime implies that anelas- 
ticity plays an essential role in the velocity-strengthening 
behavior seen in the experiments. 

In contrast, in the low / regime, the numerical simula- 
tion no longer agrees with the experiments (as shown in 
FIG.m) This indicates that a physical process comes into 
play that is not modelled in the simulations. It is known 
that the negative shear-rate dependence is due to the 
healing of the contact between grains 25|, |26| . Note that 
this process is not implemented in the numerical simula- 
tions, and this must be the reason for the disagreement. 
Thus, the frictional healing process must be properly in- 



corporated into numerical models that involve the low / 
regime. 

Taking the above points into account, we propose a 
constitutive law that holds for both the high and the low 
/ regimes. 



= cil - aln7T -I- fiQ, 



(4) 



where r is the characteristic time for frictional healing. 
The first term on the right-hand side originates from the 
dissipation due to anelasticity, and the remaining terms 
are due to the dissipation caused by intergranular fric- 
tion. The linear combination of these two effects may be 
justified if the intergranular friction force and the viscous 
force (due to anelasticity) are orthogonal. From Eq. (jH), 
the crossover inertial number Ic is given by 



a 

Cl 



(5) 



Constitutive law (U) and the critical inertial number (O 
explain the experimental data well. We conclude that 
the switch of the dominant physical processes at / ~ /c 
leads to the crossover from negative to positive shear-rate 
dependence of the friction coefficient; the anelasticity is 
dominant over the healing for I > Ic, whereas the fric- 
tional healing is essential for I < Ic- As the crossover 
is located at / ~ 10^'^, it is not surprising that it has 
not been identified in previous exp eriments on either 
inclined-plane flow (/ > 10~^) [23[ or simulated fault 
gouge (/ < 10-2) [li,[23-[33. Constitutive law (g]) and 
the resultant crossover constitute the main conclusion of 
this paper. In the following, we discuss four important 
points that are peripherally related to this main conclu- 
sion. 

First, we remark that the crossover cannot be ex- 
plained within the framework of a phenomenological con- 
stitutive law known in earthquake physics [l8|. It is 
known that a similar crossover may occur if the shear 
rate exceeds the characteristic rate of frictional healing, 
1/t 27, 2^; in this case, the steady-state friction coef- 
ficient is written as fj, — + f3liiV with /? > 0. If we 
were to adopt this relation to explain the present data 
in the high / regime, the parameter /3 would be much 
larger (0.06) than the typical values (0.005 - 0.015) that 
are reported for rocks and soda-lime glass [2^, ■ 

Second, we remark that Eq. ([3]) is different from the 
well-known relation ^ = fi^, + dH/dx, where /j,* is a con- 
stant and X is the shear displacement of a boundary. The 
latter states that the increase in the friction coefficient 
is due to the additional work required for dilation jsij. 
However, it does not state anything about the steady- 
state friction coefficient because dH/dx = at steady 
states (except for fluctuations) , whereas Eq. ^ involves 
the difference between two steady states. 

The third point we wish to discuss is with regard to 
the niechano-chemical reactions caused by frictional heat, 
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which generally reduce friction to a considerable degree 
[s^ l- We believe that anelasticity- induced strengthen- 
ing is also relevant to such systems because mechano- 
che mical rcactionij affect the nature of intcrgranular hie — 
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tioni which ia orthogonal to and independent of anclaa — 
ticity. A more detailed study on this issue is in progress. 

The fourth point of interest is the application of the 
present result to other flow geometries. Among them, 
here we derive the exponential velocity profile that is 
ubiquitously seen in heap flow experiments 12, T^]- We 
begin with the force balance equation in the flow direc- 
tion. 
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da 
dh 



pg sin 0, 



(6) 



where p is the mass density, g is the gravitational ac- 
celeration, 9 is the angle between the flow direction and 
the horizontal plane, a is the shear stress, and h is the 
depth. The existence of the body force in the flow direc- 
tion distinguishes the heap flow from the annular channel 
flow. In deriving the exponential velocity profile, it is es- 
sent ial to assume Janssen's law; i.e., the pressure within 



[14] Without such precautions, various mechano-chemical ef- 
fects (or even melting) may come into play and may 
change the nature of the granular matter 32]. 
[15] Because sliding velocity varies with radial position on the 
upper plate, we introduced an equivalent sliding velocity 
Ve = (0/3) X (Df ~Dl)/{D^-Df), where d is the angu- 
lar velocity, and Di and D2 are the inner and the outer 
diameters, respectively. With this definition and the as- 
sumption of the uniform shear stress r over the sliding 
area S, the rate of the total frictional work equals tV^S. 
The equivalent sliding velocity is referred to simply as 
the sliding velocity V throughout this study. 
The gaps between the upper plate and the side wall are 
smaller than 50 /im, which is four times smaller than the 
grain diameter. 



[16; 
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the heap, P, is independent of the depth due to the fric- 
tional force of a container. Then Eq. ^ together with 
Janssen's law leads to 



dj dp, pg sin ( 
dh dj 



P 



(7) 
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Inserting Eq. ([4]) into Eq. ((T)), we obtain 7 ~ 796" 
where A = pgsmO/aP. If Janssen's law does not hold, 
the friction coefficient is tan6' and independent of the 
dep th. Then Eq. dH) leads to a linear velocity profile. 



which is also observed in some experiments [33, 34] 
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